Anomalous dispersion of optical phonons at the neutral-ionic transition: 
Evidence from diffuse X-ray scattering 
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Diffuse X-ray data for mixed stack organic charge-transfer crystals approaching the neutral-ionic 
phase transition can be quantitatively explained as due to the softening of the optical phonon branch. 
The interpretation is fully consistent with vibrational spectra, and underlines the importance of 
electron-phonon coupling in low-dimensional systems with delocalized electrons. 



The physics of low-dimensional materials is governed 
by the complex interplay between electron-electron (e-e) 
and electron-phonon (e-ph) interactions. The resulting 
phenomenology is very rich and includes complex phase 
diagrams, competing phases, quantum phase transitions 
and multistability. Solitons and domain boundaries are 
intriguing low-energy excitations in 1-dimensional (ID) 
systems such as polyacetylene [l[ or organic charge- 
transfer (CT) salts [2J. The neutral -ionic phase transi- 
tion (NIT) observed in CT crystals with mixed stacks 
3] offers unique opportunities for investigating this rich 
physics. NIT is an electronic quantum phase transition, 
driven by the volume compression on cooling or under 
pressure, in which a stack of largely neutral molecules 
with regular spacing ...DP+ Ap~DP+ Ap~ ... (p < 0.5) 
switches to a dimerized stack of more ionic molecules, 
...(Dp+Ap-)(Dp+Ap~).... or ...(Ap-~Dp+){Ap-~Dp+) (p > 
0.5). Dimerization is a structural or Peierls transition 
driven by e-ph interactions. NIT with discontinuous 
changes in p and in the dimerization amplitude (5) are 
observed in some CT crystals, while others show contin- 
uous variations of p and 5 

In addition to domain walls or topological solitons [l[ 
that separate regions with opposite dimerization, CT 
crystals with mixed stacks have neutral-ionic domain 
walls (NIDW) that separate neutral and ionic regions 
along the ID stack. Nagaosa Q introduced the concept 
of NIDW and the corresponding excitations or domains, 
called lattice-relaxed exciton strings, LR-CT. These con- 
cepts were subsequently used to discuss the anomalies ob- 
served at NITs including the appearance of unassigned 
bands in infrared (IR) spectra [J], orders-of-magnitude 
enhancements of dielectric constants [H, 0], and diffuse 
X-ray (DXR) signals However, NIDWs were orig- 

inally restricted to the immediate vicinity of a dis- 
continuous NIT, and the energy of domains is in most 
cases too large to allow for appreciable thermal popula- 
tion Here we demonstrate that soft modes and e-ph 



coupling provide alternative explanations of all anomalies 
observed at NIT, without the restriction to a discontinu- 
ous NIT and without invoking either NIDWs or LR-CTs. 

The first observation of a soft mode associated with 
NIT was in the mid-IR spectra of tetrathiafulvalene- 
chloranil (TTF-CA) 0. TTF-CA is the prototypical 
and best characterized CT salt, with a discontinuous NIT 
at T c = 81 K from a regular neutral (p ~ 0.3) stack to 
a dimerized ionic (p ~ 0.6) stack. For T > T c weak IR 
bands polarized along the stack axis appear close to the 
frequency of Raman-active A g molecular modes. These 
bands, originally assigned to the NIDWs dynamics 0], 
are combinations of the A g modes with the IR-active 
dimerization phonon. A careful analysis of IR and Ra- 
man spectra allowed for a reliable estimate of the fre- 
quency of the soft mode [13]. Similar results have been 
obtained recently for dimethyl-TTF-CA (DMTTF-CA) 
[lH . in which almost complete softening is observed. The 
circles in Fig. Q] show the frequency of the soft mode in 
both crystals for T > T c . 

The dielectric peak at NIT was assigned to NIDW fluc- 
tuations 0, IH , but its magnitude can be quantitatively 
modeled as due to the soft dimerization mode that ac- 
quires a huge IR intensity near the transition 12, [l3j ]. 
DXR data @, |j| have been interpreted in terms of long 
LR-CTs in DMTTF-CA and short ones in TTF-CA: so 
far they offer the only unchallenged evidence for NIDWs. 

DXR scattering was widely used in the 80's to study 
soft modes that drive structural transitions, includin g th e 
Peierls transition in CT salts with a segregated stack [lj| . 
The Peierls transition in mixed (DA) stacks is however 
different: the soft mode or incipient soft mode is an op- 
tical phonon instead of an acoustic phonon as in segre- 
gated stacks. The sharp Lorentzian-shaped DXR signal 
collected in one spatial direction implies a strong disper- 
sion of the relevant phonons and is unusual for the opti- 
cal branch of a molecular crystal. This was the basis for 
rejecting a soft mode interpretation in favor of LR-CTs 
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FIG. 1: Circles: the soft mode frequency from refs. [lOl. [Tlj. 
Diamonds: the soft mode frequency (arbitrary units) as in- 
ferred on the basis of Eq. [T]from DXR data in Ref. @, H|- 



[8]. But e-ph coupling is extremely effective in ID and 
anomalous dispersion of optical phonons is well known in 
polyacetylene [15] . 

Here we demonstrate that DXR data can be quantita- 
tively understood as due to the progressive softening of 
the zone-center optical phonons as the system is driven 
towards the NIT. A sharp anomaly develops in the optical 
phonon branch that, although occurring in non-metallic 
systems, is analogous to the Kohn anomaly in ID met- 
als, and shares the same origin, namely, e-ph coupling. 
Peierls transitions were similarly proposed for ID metals 
and later generalized to transitions driven by e-ph cou- 
pling in insulators and spin chains. 

The DXR intensity, /, is related to the frequency u of 
the soft mode by a simple relation [lij ]: 
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where q is the wavevector. Reliable estimates of w(0) 
were obtained from the analysis of combination bands in 
IR spectra of both TTF-CA and DMTTF-CA [n], El. 
Since DXR scattering is in arbitrary units, Eq. [1] yields 
the temperature dependence of w(0) to within a scale 
factor. The diamonds in Fig. Q] are scaled w(Q) from 
DXR. The striking coincidence of DXR and IR estimates 
of w(0) for T > T c in both salts provides strong moti- 
vation for a soft-mode analysis. But several issues must 
be addressed for such an interpretation, starting with the 
strong dispersion of the optical branch. 

To confirm the soft-mode picture, we consider the dis- 
persion of optical phonons in a ID stack with equal spac- 
ing, harmonic force constant K for nearest-neighbor lat- 
tice displacements, molecular masses for D and A, and 
linear e-ph coupling. For simplicity, we describe the 
electronic structure of the mixed stack using a spinless- 
fermion (SF) model, an uncorrelated model [9( that can 
readily be solved for long chains. The electronic Hamil- 
tonian describes N/2 spinless fermions on N molecular 
sites with periodic boundary conditions, as follows: 



H, 



SF 



= r£(- 



i) 



(2) 





-| go v 






75 




(cm 


50. 




3 








25 






0> 
1 






r 


0.5 




FIG. 2: The T-dependence of the optical dispersion curve 
to(q), calculated for parameters relevant to DMTTF-CA, e<i = 
0.4 and loo = HO cm^ 1 on 902 sites. The continuous red line 
is the dispersion curve for = 0. 



where p runs on the sites, h p counts the fermions on p 
site, t p is a transfer integral, and b p = (a p a p +i + H.c.) 
is the bond-order operator. In the SF model, a single 
fermion moves from an odd (D) to an even (A) site to 
describe electron transfer between DA and D + A~. 2r 
then represents the energy required to ionize a DA pair 
and implicitly includes all contributions from e-e interac- 
tions and from the coupling of electrons with molecular 
vibrations Q. The hopping integrals are modulated by 
lattice phonons: t p = tQ+a(x p —x p+ i) where x p is the dis- 
placement of the p site, and t will be used as the energy 
unit. The e-ph coupling constant a enters the definition 
of the vibrational relaxation energy: td = a 2 /K . Aside 
from a sp in factor, the model was discussed by Rice and 
Mele [Hf, and at V = it reduces to the SSH model Q. 

The vibrational problem in the absence of e-ph cou- 
pling (e^ = 0) describes a ID harmonic chain with two 
masses. It has a well known solution, and the red con- 
tinuous line in Fig. [2] shows the dispersion curve for the 
optical branch, u>o(<7)j calculated for the molecular masses 
of DMTTF-CA and with K fixed as to get w (0) = 110 
cm . In the presence of e-ph coupling, td > 0, new 
harmonic force constants add to K, corresponding to the 
second derivatives of the ground state (GS) electronic 
energy E vs the site displacements: 
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where n r = — (d 2 E /dt p dt p+r ) eq , the bond-bond polariz- 
abilities, are electronic quantities. In the SF model: 
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where \kl) is the excited state of Hsf with a fermion 
promoted from the fc-th filled orbital to the l-th empty 
orbital, and ek is the fc-orbital energy. Much as in 
polyacetylene 17] . the electronic derealization leads to 
long-range bond-bond polarizabilities, and hence to long- 
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range force constants that are responsible for an anoma- 
lous dispersion of the optical phonon branch. The acous- 
tic branch is only marginally affected by e-ph coupling. 

In momentum space, the squared vibrational frequen- 
cies Lo(q) 2 are obtained, as usual, by diagonalizing the 
2x2 force constant matrices, whose elements contain 
sums over II r as well as K and the molecular masses. 
The largest softening occurs at q = 0, where a simple 
expression holds for the optical branch frequency [3]: 



DMTTF-CA 

6000 •, 



Lj(0) = Lj o (p)y/l-e dX (r) 



(5) 



where x — — {d 2 E(T, S)/d6 2 )o, the curvature of the elec- 
tronic GS energy along the dimerization coordinate, S, 
only depends on T. The divergence of \ a t T = marks 
the unconditional Peierls instability of a half-filled ID 
metal. For finite €d CT stacks are stable against dimer- 
ization for r > Tp, where Tp is defined by Cdxi^p) = 1- 
Peierls transitions have been extensively studied and 
there are analytical expressions for for uncorrclatcd 
models and numerical results for correlated models 
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Figure [2] shows the evolution of u>(q) as function of 
r > Tp. For the chosen parameters, the Peierls transi- 
tion occurs at Tp = 0.05, where the softening is complete, 
u>(0) = 0, and a giant anomaly develops in the dispersion 
of optical phonons. The softening of the optical branch 
follows directly from linear e-ph coupling. An important 
new result in Fig. [5] is the explicit expression for the 
phonon dispersion as the dimerization transition is ap- 
proached. With w(q) in hand, we can model the DXR 
profile 7(g) using Eq. [TJ 

DXR profiles I(q) are measured as a function of T, 
while the phonons in the SF model are computed as a 
function of T. The profile in either case is characterized 
by the peak height 1(0) and by the half-width at half- 
maximum, qi/2- The calculated I(q) in Fig. [3] are based 
on Eq. [H with lu(0) taken from vibrational data, and the 
best- fit parameters, e<j and ljq(0), in the caption. Insets 
b and e specify T(T) for DMTTF-CA and TTF-CA, re- 
spectively, and an almost linear relation is extracted in 
each case. The measured and calculated values of qi/ 2 
are seen in insets c and / to be in excellent agreement, 
especially close to the transition. DXR profiles in mixed- 
stack CT salts close to the NIT can be quantitatively 
explained by the softening of the dimerization mode and 
the related anomaly in the optical phonon branch. 

The sharp DXR peaks observed in DMTTF-CA were 
previously assigned to the presence of long LR-CT exci- 
tations [2], whereas the broader signals in TTF-CA were 
ascribed to shorter domains Q. Here the difference be- 
tween the two systems is quite naturally related to their 
different transitions. DMTTF-CA undergoes a continu- 
ous (or almost continuous) dimerization transition 111 ]: 
p(T) is continuous, the phonon frequency ui(0) softens 
to zero (or almost so), and the dip in the phonon dis- 
persion fully develops leading to sharp DXR peaks. On 
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FIG. 3: Main panels (a, d): the T-evolution of the DXR 
signal as calculated in the SF model. Panel a: DMTTF-CA, 
e d = 0.4, oj (0) = 110 cm" 1 ; panel d: TTF-CA, e d = 0.55, 
wo(0) = 130 cm -1 . Insets 6 and e show the T(T) relation; c 
and / compare the width of calculated DXR peaks (circles) 
with the experimental values (crosses) from refs @, @|- 



the other hand, only incipient soft-mode behavior is ob- 
served in TTF-CA, whose discontinuous NIT at 81 K 
interrupts the softening As a result, the dispersion 
anomaly does not fully develop and comparatively broad 
DXR peaks are observed. 

The non-interacting model, H$f in Eq. [21 can be 
solved for very long chains, as needed to get enough 
resolution in the dispersion curves. To show that the 
anomalous dispersion in the optical phonon branch sur- 
vives correlations, we now turn attention to the standard 
correlated model for NIT, the modified Hubbard (MH) 
model Q with linear e-ph coupling. MH describes N 
electrons on N sites, and, as in SF, 2r is the energy re- 
quired to ionize a DA pair. Correlations are taken into 
account by a restricted basis that excludes D 2+ and A 2 ~ 
sites. Intersite e-e interactions V can be introduced in 
addition to linear e-ph coupling. 

MH bond-bond polarizabilities n r are calculated as nu- 
merical derivatives of the exact GS energy of finite chains 
(up to 20 sites) with periodic boundary conditions. We 
compare results for SF and MH chains with the same p, 
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FIG. 4: Left panels: r-dependence of the bond-bond polariz- 
abilities, calculated for p=0.46 (dashed line for SF, triangles 
for MH) and 0.22 (continuous line for SF, circles for MH). 
Right panel: the dispersion curves calculated for the same p 
as in the left panels, and parameters relevant to DMTTF- 
CA. Lines and symbols refer to SF (N = 902) and MH model 
{N = 20). 



often been interpreted as evidences for exotic excitations 
like NIDW or LR-CT. Such excitations are, however, only 
expected in the immediate vicinity of a first-order tran- 
sition 0, Q ■ The soft mode interpretation of anomalous 
mid-IR bands in CT salts appearing on approaching NIT, 
as well as the huge peak in the dielectric constant mea- 
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sured at NIT has already been discussed 
Here we showed that DXR scattering is another manifes- 
tation of e-ph coupling in ID. 

Work in Parma supported by NE MAGMANET 
NMP3-CT2005-515767 and Italian MIUR through FIRB- 
RBNE01P4JF. 



and in the left panels in Fig. 2]we show the r-dependcnce 
of the n r for two systems with p = 0.22 and 0.46, far and 
close to NIT, respectively. For both systems, apart from 
a scale factor, the 11^ calculated in SF and MH models 
have similar behavior and oscillate between positive and 
negative values, much as found in polyacetylene 17[. The 
right panel of Fig. 0] compares the MH and SF phonon 
dispersion curves calculated for the same p values. For 
the SF case we set the parameters relevant for DMTTF- 
CA. In the MH model, to account for the stronger e-ph 
coupling in correlated systems, we rescale ed to obtain 
the same w(0) as in SF. 

Far from NIT, at p — 0.22, the effect of e-ph coupling 
is mild: n r decays rapidly with r, well within distances 
achievable in a periodic 20 site chain. As a result, the 
discrete MH points for uj(q) fall exactly on top of the SF 
curve in Fig. 4, right panel. On the contrary, near the 
transition (p = 0.46), n r acquires a long-range tail that 
cannot be appreciated in a 20-site system. Nevertheless, 
the MH points of the dispersion curves approximately fall 
on the SF curve. We conclude that the sharp anomaly 
in the optical phonon branch is a characteristic feature 
of NIT, and that its basic physics is well captured by 
models without e-e correlation. 

In conclusion, we ascribe DXR scattering in DMTTF- 
CA and TTF-CA crystals on approaching the NIT to 
the evolution of a Kohn-like anomaly in the optical 
phonon branch, related to the softening of the dimeriza- 
tion phonon. The anomaly is due to amplified e-ph cou- 
pling near an electronic transition. Our interpretation, 
quantitatively consistent with vibrational frequencies as 
extracted from mid-IR data, is supported by the calcula- 
tion of phonon dispersion curves in SF and MH models 
for electronic systems with linear (Peierls) e-ph coupling. 
Linear e-ph coupling has enormous effects in systems 
with delocalized electrons in ID, as best demonstrated 
by polyacetylene [3, [l|| . In CT salts the amount of der- 
ealization increases when moving towards NIT from the 
N side, giving rise to several anomalous features that have 
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